Regularity and chaos at critical points of first-order quantum phase transitions 
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We study the interplay between regular and chaotic dynamics at the critical point of a generic first 
order quantum phase transition in an interacting boson model of nuclei. A classical analysis reveals 
a distinct behavior of the coexisting phases in a broad energy range. The dynamics is completely 
regular in the deformed phase and, simultaneously, strongly chaotic in the spherical phase. A 
quantum analysis of the spectra separates the regular states from the irregular ones, assigns them 
to particular phases and discloses persisting regular rotational bands in the deformed region. 
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Quantum phase transitions (QPTs) at zero tempera- 
ture arc structural changes of the system's ground state 
resulting from a variation of parameters A in the quan- 
tum Hamiltonian H{X) Aspects of QPTs and their 
finite- A'' precursors are currently attracting considerable 
theoretical and experimental interest, as they occur in 
diverse dynamical systems including spin lattices J2| , en- 
sembles of ultracold atoms [3j] and atomic nuclei [J]. 

The particular type of QPT is reflected in the topology 
of the underlying mean-field (Landau) potential ^(A). In 
first order (discontinuous) QPTs, V{X) develops multiple 
minima, corresponding to distinct phases (characterized 
by distinct values of an order parameter) that coexist in 
a range of A-values and become degenerate at the crit- 
ical point A — Ac. In higher order (continuous) QPTs, 
the single minimum of ^(A) shifts continuously but non- 
analytically to a different value of the order parameter at 
A = Ac, without phase coexistence. 

Interacting many-body systems undergoing QPTs pro- 
vide a fertile ground for studying the emergence of quan- 
tum chaos. A typical example, is the recent analysis in 
quantum optics models of A^ two-level atoms interact- 
ing with a single-mode radiation field [1], where the on- 
set of chaos is triggered by continuous QPTs. In this 
communication, we turn our attention to systems with 
discontinuous QPTs, where multiple- well potentials and 
phase coexistence provide even a richer environment to 
study the interplay of order and chaos. Indeed, studies 
of 2D single-particle models indicate that the classical 
motion can be dissimilar (regular or chaotic) in differ- 
ent local minima f^. In the present work, we show that 
such mixed form of dynamics and non-uniform onset of 
chaos arises from rotational-invariant Hamiltonians in a 
system of A' interacting constituents, undergoing a first- 
order QPT. A marked separation between regular and 
chaotic dynamics is observed in both classical and quan- 
tum analyses. 

As a concrete example, we adopt the interacting bo- 
son model (IBM) Q, widely used in the description of 
quadrupole collective states in nuclei, in terms of a sys- 
tem of A monopole (s) and quadrupole (d) bosons, repre- 
senting valence nucleon pairs. The total boson number N 
and angular momentum L are conserved, which enables 
to solve the model exactly by numerical diagonalization 



of finite matrices for any fixed N and L. The model 
accommodates both continuous and discontinuous quan- 
tum shape-phase transitions [8], which are manifested 
empirically between spherical and deformed nuclei [4|. 

Apart from kinetic rotational terms, the most general 
two-body IBM Hamiltonian at the critical point of a first 
order QPT, can be transcribed in the form Q 
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P2(/3o)/A(A-l) 

(2) 
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where PliM = l3osUl + ^7/2 {d^d%'\ P^M = 
(— l)^P2.-/i(/3o) and the dot implies a scalar product. To 
facilitate the comparison with the classical limit, involv- 
ing large A^, the critical-point Hamiltonian ([IJ is scaled 
by A(A^ — 1). For f3o > 0, Hch annihilates two distinct 
degenerate ground states, which have the form of static 
condensates \I3;N) = {m)-^/^Tl{l3)^\0) of A^ bosons, 
rt(/3) = [pdl + V2 - /32st]/V2, with (i) 13^0 (spherical 
s-boson condensate) and (n) /3 = V2/3o(l -I- /3o)~"^^^ (de- 
formed condensate) and correspond to the two coexisting 
shape phases of the nucleus. The parameter /3o in Hen 
determines the equilibrium deformation in the deformed 
phase as well as the height of the potential barrier be- 
tween the two phases (see below). 

The classical limit of the IBM is obtained through the 
use of coherent states [l3|. This amounts to replacing 
(fit, dj^) by six c- numbers (a*, a*) rescaled by ^/N and 
taking N ^ oo, with 1/A^ playing the role of h. Number 
conservation ensures that phase space is 10-dimensional 
and can be phrased in terms of two shape (deformation) 
variables, three orientation (Euler) angles and their con- 
jugate momenta. Chaotic properties of the IBM have 
been studied extensively [ll|, albeit, with a simplified 
Hamiltonian, giving rise to an extremely low, hence non- 
generic, barrier. The Hamiltonian of Eq. ([1]) can accom- 
modate a barrier with an adjustable height, hence allows 
a systematic study of its impact on the dynamics. In the 
classical analysis presented below we consider, for sim- 
plicity, the dynamics of L = vibrations, for which only 
two degrees of freedom are active. The rotational dynam- 
ics with L > is examined in the subsequent quantum 
analysis. 

For the particular case of the Hamiltonian ([T]), con- 
strained to L = 0, the above procedure yields the follow- 
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ing classical Hamiltonian 

+/3o^(l-Hd,o)/2 
X [(p^//3 ~ (3pl - (3^) cos 37 + 2ppp^ sin 87] . (2) 

Here the coordinates /3 £ [0, \/2], 7 € [0, 27r) and their 
canonically conjugate momenta Pjs € [0, \/2] and pj G 
[0, 1] span a compact classical phase space. The term 
Tidfi = (r+^2)/2, with T ^ pj+p?^/l3^, denotes the clas- 
sical limit of the d-boson number operator fid — '^jj'^A' 
(restricted to i = 0) and forms an isotropic harmonic os- 
cillator Hamiltonian in the (3, 7 variables. Notice that ([2]) 
contains complicated momentum-dependent terms origi- 
nating from the two-body interactions in ([1}, not just the 
usual quadratic kinetic energy T. Setting = = in 
Eq. ([2|) leads to the potential 

Vcri/h2 = i/32/?2 + i(l-/32)/34 

-i/?ov/2-/32/33cos37 , (3) 

which can be alternatively obtained for 7 = as an 
expectation value of -ffcri © in the static condensate 
mentioned above, Vcri(/3,7 = 0) = {l3; N\Hcri\l3; N) . 
Apart from being compact, the /3,7 variables are anal- 
ogous to the coordinates of the geometric collective 
model [HI and portray the amount of quadrupole de- 
formation (j3 — corresponds to the spherical, while 
/3 > to deformed shapes) and the triaxiality of the 
nucleus (7 — 27rfc/3, with fc = 0, 1,2 correspond to pro- 
late and 7 = tt/3 + 27rfc/3 to oblate axially-symmetric 
shapes, while other values to triaxial shapes). They may 
be interpreted as polar coordinates in an abstract plane 
parametrized by Cartesian coordinates x = P cos 7 and 
y = /3sin7. The potential of Eq. ([3|) displays a three- 
fold rotational symmetry about the origin /3 = 0. Its 
four minima, being degenerate at energy Knin — 0, are of 
two types and correspond to the different phases: (i) the 
single spherical minimum ai (3 — and (ii) the three 
equivalent deformed minima at /3 = -\/2/3o(l + /^o)"^^^ 
with 7 = 0, 27r/3, 47r/3. In between the spherical and de- 
formed minima, there are three equivalent saddle points 
located at /? = [1 - (1 + /Sq)"^^^]^^^' which create po- 
tential barriers of height Vb = /i2[l - (1 + Plf'^fl'^, 
separating the two phases. The limiting value at the do- 
main boundary is V[(3 = \/2, 7) = /12. The potential of 
Eq. dHl) is similar to that considered in 0], but differs by 
the square-root term and the compact domain. 

The nature of the classical motion associated with 
the Hamiltonian ([2]) can be depicted conveniently via 
Poincare surfaces of sections in the plane y = 0, plotting 
the values of x and the momentum px each time a trajec- 
tory intersects the plane. Regular trajectories are bound 
to toroidal manifolds within the phase space and their 
intersections with the plane of section lie on ID curves 
(ovals). In contrast, chaotic trajectories randomly cover 
kinematically accessible areas of the section [l^ . 



Poincare sections for (3q = 1.0,1.3,1.5, are displayed, 
respectively, in the left, middle and right columns of 
Fig. [H The three cases correspond to "low" , "medium" 
and "high" potential barriers Vb/h2 = 0.04, 0.10, 0.16, 
(compared to Vfc//i2 — 0.0009 in previous works [11]). In 
each case (see panels ak-bk-Ck), we show the sections at 
five different energies Ek (k =1-5) relative to the barrier- 
tops Vb-. El = 0.2514, E2 = 0.9Vb are below, E3 = Vb 
at, and E4 = l.lVb, E^ = 4Vb are above the barrier. 
The bottom row (panels a-b-c) depicts the correspond- 
ing classical potential Vc,:i{x, y — 0), Eq. ([3]). In all three 
cases, the motion is predominantly regular at low energies 
and gradually turning chaotic as the energy increases. 
However, the classical dynamics evolves differently in the 
vicinity of the two wells. The family of regular trajecto- 
ries near the deformed minimum has a particularly simple 
structure. It forms a single set of concentric loops around 
a single stable (elliptic) fixed point. The trajectories re- 
main regular even at energies far exceeding the barrier 
height Vb- In contrast, a more complicated structure de- 
velops near the spherical minimum. Here, at low energies 
(panels ai-61-ci), one observes four major islands sur- 
rounding stable fixed points, and unstable (hyperbolic) 
fixed points in-between. As the energy increases (approx- 
imately for i? > Vfc/3 [l3|), considerable fraction of the 
trajectories becomes chaotic (panels 02-62-02) until com- 
plete chaoticity is reached near the barrier top. The clear 
separation between regular and chaotic dynamics, associ- 
ated with the two minima, persists all the way to i? = 14, 
where the two regions just touch (panels 03-63-03). At 
E > Vb, the chaotic trajectories from the spherical re- 
gion can penetrate into the deformed region and a layer 
of chaos develops (panels 04-64-04), and gradually dom- 
inates the surviving regular island for i? ^ Vf, (panels 
05-65-05). In general, the regularity is more pronounced 
for higher barriers (larger /3o). 

Quantum manifestations of classical chaos are often 
detected by statistical analyses of energy spectra [T^. 
In a quantum system with mixed regular and irregular 
states, the statistical properties of the spectrum are usu- 
ally intermediate between the Poisson and the Gaussian 
orthogonal ensemble (GOE) statistics. Such global mea- 
sures of quantum chaos are, however, insufficient to re- 
flect the rich dynamics of an inhomogeneous phase space 
structure, as in Fig. 1, with well separated regular and 
chaotic regions. To do so, one needs to distinguish be- 
tween regular and irregular subsets of eigenstates in the 
same energy intervals. For that purpose we employ the 
spectral lattice method of Peres [T^ , which provides ad- 
ditional properties of individual energy eigenstates. The 
Peres lattices are constructed by plotting the expectation 
values Oi = {i\0\i) of an arbitrary operator, [0,H] ^ 0, 
versus the energy Ei = {i\H\i) of the eigenstates The 
lattices {Oi,Ei} corresponding to regular dynamics can 
be shown to display a regular pattern, while chaotic dy- 
namics leads to disordered meshes of points. The method 
has been recently applied to quantum chaos in the geo- 
metric collective model of nuclei [16| . 
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FIG. 1: (Color online). Poincare sections for the classical 
Hamiltonian ((2]) with /12 = 1 and /3o = 1.0, 1.3, 1.5 (left, mid- 
dle and right columns, respectively). In each case, the surfaces 
are drawn at energies Ek {k = 1-5) relative to the barrier- 
tops Vb-. El = 0.25V6, E2 = 0.9V6, E3 = Vi, E4 = l.lVi,, 
E5 — 4Vb (panels ak-b^-Ck). The bottom row (panels a-b-c) 
depicts the corresponding classical potential Vcri{x,y — 0), 
Eq. ([Sjl, with coexisting spherical and deformed minima at 
a; = and i > 1, respectively. The reference energies, E^, 
are indicated by horizontal dashed and dotted lines. Vertical 
lines mark the turning-points for E^ < 14. Notice the pro- 
nounced separation of regular and chaotic regions in phase 
space, associated with the two minima. 



In the present analysis we choose the Peres operator to 
he O = fid- The lattices correspond to the set of points 
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FIG. 2: (Color online). Quantum Peres lattices {xi,Ei} for 
L — eigenstates \i) of _ffcri (O, with h2 — 1, N — 80 and 
Po = 1.0, 1.3, 1.5. The quantity Xi — ^y2{i\'hd\i) /N enables a 
direct comparison with the classical potential Vcii{x,y = 0), 
Eq. ([3}. The reference energies Ei, E2, E4, Es, are indicated 
as in Fig. [1] Notice the sequences of regular states in the 
deformed region. The lowest sequence consists of L = 
bandhead states of the ground g{K — 0) and I3"{K = 0) 
bands. Regular sequences at higher energy correspond to 
P"-f'^{K = 0), /3"7''(_K' = 0) bands, etc. 



{xi, Ei}, with Xi = y^2{i\nj^\i) /N and \i) being the eigen- 
states of Hcii (IIJ . The expectation value of hd in the con- 
densate, X =^ l3 = y/2{l3;N\nd\(3;N)/N, is related to the 
deformation /3 (whose equilibrium value is the order pa- 
rameter of the QPT) and the coordinate x in the classical 
potential, Vcriix,y = 0) = Vcd(/3,7 = 0), Eq. ©. The 
spherical ground state has rid = Xi =0. Excited spherical 
states show characteristi c domi nance of single compo- 
nents 0, hence Xi ^ ^/ud/N is small for Ud/N « 1. 
Rotational members of the deformed ground band are 
obtained by L-projection from |/3; N) and have Xi ~ f5 
to leading order in N. This relation is still valid, to 
a good approximation, for states in excited deformed 
bands, whose intrinsic states are obtained by replacing 
condensate bosons in |/3; N) with orthogonal bosons rep- 
resenting 13 and 7 excitations 17]. These attributes have 
the virtue that the lattices {xi,Ei} can identify the reg- 
ular/irregular quantum states and associate them with a 
given region in the classical phase space. 

Fig. [5] presents the lattices calculated for L ^ Q eigen- 
states of Hcii with N — 80 and /3o — 1-0, 1.3, 1.5. In each 
case, one can clearly identify regular sequences of states 
localized within and above the respective deformed wells. 
They form several chains of lattice points close in energy, 
with the lowest chain originating at the deformed ground 
state. A close inspection reveals that the x^-values of 
these regular states, lie in the intervals of values oc- 
cupied by the regular tori in the Poincare sections in 
Fig. [TJ Similarly to the classical tori, these regular se- 
quences persist to energies well above the barriers Vt- In 
contrast, the remaining states, including those residing in 
the spherical minimum, do not show any obvious patterns 
and lead to disordered (chaotic) meshes of points at high 
energy E > Vb- (Although there are only a few lattice 
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FIG. 3: (Color online). Quantum Peres lattices {xi,Ei} for 
eigenstates of Hcvi (HI with h2 = 1, N — 50, /3o = 1.5 and 
L = 0, 2, 3, 4. To enhance visibility, a small term is 
added, leading to an energy shift AEl = 0.0032L(L + 1), 
without affecting the wave functions. The classical poten- 
tial Vcri{x,y = 0) and reference energies Ek are indicated as 
in Fig. [21 Notice the rotational bands {K = 0, L = 0, 2, 4) 
and {K = 2, L = 2, 3, 4) formed by the regular states in the 
deformed phase. 

points near x = 0, the chaotic behavior of these spherical 
states is evident by glancing at the Poincare sections of 
Fig. [J). The ability of the Peres method to separate regu- 
lar states from irregular ones can be tested by fitting the 
nearest neighbors level spacing distribution of the quan- 
tum spectrum by a Brody distribution Puj{S) [l^. The 
parameter ui interpolates between the Poisson (w = 0) 
and GOE statistics (a; = 1), corresponding to integrable 
and fully chaotic classical motion, respectively. We deter- 
mine uj for the full (mixed) spectrum of L = eigenstates 
and for a partial ("chaotic") spectrum, from which the 
tentatively regular levels are excluded. Specifically, for 
{N = 80, Po = 1.5), we eliminate the 32 levels which 
form the regular sublattice in the deformed phase (see 
Fig. Ht) from a set of 364 states with E < 0.8. For the 
mixed spectrum we obtain uj = 0.78 ± 0.10 and for the 
"chaotic" spectrum ui ~ l.OliO.lO, thus confirming the 
regular character of the excluded levels. 

The Peres lattices corresponding to eigenstates of 



Hcri © with L = 0,2,3,4, TV = 50 and /3o = 1.5 are 
shown in Fig. [31 They disclose families of regular pat- 
terns in the deformed phase. These correspond to rota- 
tional bands of states having a common intrinsic struc- 
ture, as indicated by their nearly equal values of {fid). 
The regular L — states form bandheads of rotational 
sequences L = 0, 2, 4, . . . [K = bands) and are accom- 
panied by sequences L = 2,3,4,... built upon regular 
L = 2 states {K = 2 bands). Additional iiT-bands with 
L — K, K +1, K + 2, . . . (not shown in Fig. [3]) can also be 
identified. In the nuclear physics terminology, the lowest 
K = band refers to the ground band and excited K- 
bands correspond to multiple /3 and 7 vibrations about 
the deformed shape with angular momentum projection 
K along the symmetry axis. Detailed analysis reveals 
additional band-like properties of these states, in par- 
ticular, the coherent decomposition of wave functions in 
the rotor basis 14]. Such ordered band structures are 
not present in the chaotic parts of the lattice and their 
persistence in the spectrum is unexpected, in view of the 
strong mixing and abrupt structural changes taking place 
at the critical point. It validates the recently proposed 
regularity-induced adiabatic separation of intrinsic and 
collective dynamics [l^, for a subset of states. 

In summary, we have examined the interplay of regu- 
larity and chaos in an interacting system modeling the 
dynamics at the critical point of a first order QPT be- 
tween spherical and deformed nuclei. A classical analysis, 
using Poincare sections, shows that the two coexisting 
phases exhibit different susceptibilities towards the on- 
set of chaos. While the deformed phase displays robustly 
ordered and rather simple dynamics, the spherical phase 
shows strongly chaotic behavior. An analysis of the quan- 
tum spectrum, using Peres lattices, identifies unexpected 
regular rotational bands in the deformed region, which 
are absent from the disordered portions of the lattice. 
Although the regular and chaotic motions coexist in a 
broad interval of energies, they are well separated and can 
be distinguished by both classical and quantum calcula- 
tions. The system in the domain of phase coexistence, 
thus provides a clear cut demonstration of the classical- 
quantum correspondence of regular and chaotic behavior, 
illustrating Percival's conjecture concerning the distinct 
properties of regular and irregular quantum spectra [TH . 
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